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Current strategies of indirect Dark Matter detection with neutrino telescopes are 
based on the search for high-energy neutrinos from the Solar core or from the center of 

o ' 

CN ' the Earth. Here, we propose a new strategy based on the detection of neutrinos from 

Dark Matter annihilations in mini-spikes around Intermediate Mass Black Holes. 
Neutrino fluxes, in this case, depend on the annihilation cross-section of Dark Matter 
fN| ' particles, whereas solar and terrestrial fluxes are sensitive to the scattering cross- 

ed ■ section off nucleons, a circumstance that makes the proposed search complementary 

> , 

to the existing ones. We discuss the prospects for detection with upcoming under- 
water and under-ice experiments such as ANTARES and IceCube, and show that 
several, up to many, sources could be detected with both experiments. A kilometer 
scale telescope in the Mediterranean appears to be ideally suited for the proposed 
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Current strategies of indirect Dark Matter (DM) detection with neutrino telescopes are 
based on the search for high-energy neu ;rinos from the Solar core or from the center of he 
Earth (for recent reviews see Refs. Q, Q, 3]). The idea, proposed more than 20 years ago 
is quite simple: if we live immersed in a halo of DM, then there is a wind of DM particles 
passing through the Sun and the Earth. A fraction of all these particles is expected to lose 
enough energy, by scattering off the nuclei of the star or planet, to become gravitationally 
bound. Captured particles will then slowly sink at its center, due to subsequent scattering 
events, and will annihilate producing a steady flux of high energy neutrinos (see Refs. l| 
and [3j], and references therein). 

At equilibrium, the flux of neutrinos is determined by the number of DM particles cap- 
tured per unit time by the Sun or Earth, a quantity that depends on the local density (and 
velocity distribution) of DM and on the scattering cross-section of DM particles off nucle- 



ons, i.e. the same quantities that are probed by direct detection experiments (e.g. Ref. 
and references therein). Unfortunately, a large portion of the DM parameter space will not 
be accessible to neutrino telescopes in these scenarios. In particular, only the most opti- 
mistic neutralino models in the Minimal Supersymmetric Standard Model will be accessible 
to kilometer sized experiments (e.g. Ref. [6J]), while for Kaluza-Klein DM candidates, and 
more specifically for the B^, first excitation of the hypercharge gauge boson in theories 



with Universal Extra Dimensions 
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, , between 0.5 and 10 events per year are anticipated in 
kilometer scale neutrino telescopes 8||. 

We propose here a new strategy for the search of DM with neutrino telescopes, based 
on the detection of secondary neutrinos produced by DM annihilations in 'mini-spikes', i.e. 
recently proposed enhancements of DM around Intermediate Mass Black Holes (IMBHs). 
We argue that mini-spikes can be bright sources of neutrinos, and we study the p rosp ects 
for detecting them with upcoming neutrino telescopes ANTARES j^J and IceCube jl(|. 

Mini-spikes result from the reaction of DM mini-halos to the formation or growth of 
IMBHs. In fact, the adiabatic growth of a massive object at the center of a power-law 
distribution of DM with index 7, induces a redistribution of matter into a new power- 
law with index 7 sp = (9 — 27)/ (4 — 7) 11]. This formula is valid over a region of size 
R s pa 0.2tbh, where tbh is the radius of gravitational influence of the black hole, defined 
implicitly as M(< tbh) — Mbh, with M(< r) mass of the DM distribution within a sphere 
of radius r, and Mbh mass of the Black Hole jl^ . 

To make quantitative predictions, we focus on a specific IMBHs formation scenario jl^j ]. 
representative of a class of models where these objects form directly out of cold gas in early- 
forming DM halos, and are characterized by a large mass scale, of order 10 5 M Q (see also 
Ref. [14J and references therein). In this case, mini-spikes can be bright sources of DM 
annihiiation rad,ation; indirect searches with gamma-ray telescopes, ami in oart.cu.ar with 



GLAST |l5j . expected to be launched in 2007, are particularly promising [16]. However, 
GLAST will be sensitive to energies up to 300 GeV, which means that the most interesting 
feature of the annihilation spectrum, i.e. the presence of an identical cut-off at an energy 
equal to the DM particle mass, would be missed if the DM particle is heavy. Air Cherenkov 
Telescopes (ACTs) such as CANGAROO [17], HESS Q, MAGIC [19] and VERITAS 2(J 
could extend the observations to higher energies, but due to their narrow field-of-view, they 
are not optimized for full-sky surveys. 
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FIG. 1: Sky map in equatorial coordinates showing the position of Intermediate Mass Black Holes 
in one random realization (red diamonds), and in all 200 realizations (blue dots). The concentration 
at negative declinations corresponds to the position of the Galactic center (black open diamond). 

We show here that neutrino telescopes can extend full-sky searches at higher energies . 
In case of detection, the better sensitivity and energy resolution of ACTs can then be used 
to effectively discriminate mini-spikes from ordinary astrophysical sources. 

We generate here a mock catalog of IMBHs, and associated mini-spikes, from the (200) 
realizations of Milky- Way like halos presented in Ref. Q|, by randomly assigning to every 
object of given galactocentric distance r GC , 2 angular coordinates (9, 0), following a uniform 
probability distribution on the surface of the unit sphere. The specifics of the Monte Carlo 



and the method for populating black holes at 
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halo evolution procedure are given in Ref. |21 
high-redshift are described in detail in Refs. 



22j. The resulting distribution is shown in 
Fig. [H where we chose to adopt equatorial coordinates, to emphasize the different regions 
that will be accessible to Antares and Icecube. We multiply all neutrino fluxes by a visibility 
factor V^(5), which keeps into account that the effective area of a neutrino telescope at 
terrestrial latitude 0, drops abruptly for sources above the horizon. The visibility factor is 
thus proportional to the fraction of time spent by a given source of declination 5 below the 
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horizon. By means of simple geometrical considerations we find that 



0, if |<J| >tt/2- \(f>\,6-(p > 

arccos [tan(0) tan(<5)] /vr, if |5| < tt/2 - |0| (1) 

1, if \8\ >tt/2- <0 



For Antares = 42°50', while for IceCube = —90°. In the case of IceCube the visibility 
factor reduces to 

f 0, if 8 < 

^IceCube — \ (2) 

[ 1, if 5 > 

Mini-spikes would be copious sources of neutrinos, which can be produced either by the 
direct annihilation of DM to neutrino pairs, or through fragmentation and decay of secondary 
particles such as quarks, leptons and gauge bosons. If N ye (E) is the spectrum of neutrinos 
z/£, of flavor £ — e,fi, r, per annihilation, the flux of neutrinos from an individual mini-spike 
in absence of oscillations can be expressed as [16] 

= 00^-^(^)26^^^ (E) (3) 

with 0o = 9 x 10 _10 cm _2 s _1 . The first two factors depend on the particle physics parameters, 
viz. the mass of the DM particle in units of 100 GeV m x ioo, and its annihilation cross 
section in units of 10 _26 cm 3 /s, (o~v)26, while the third factor accounts for the flux dilution 
with the square of the IMBH distance to the Earth in kpc, -Dkpc- Finally, the normalization 
of the flux is fixed by an adimensional luminosity factor L sp , that depends on the specific 
properties of individual spikes. In the case where the DM profile before the formation of the 

, the final DM 



IMBH follows the commonly adopted Navarro, Frenk and White profile 
density p(r) around the IMBH will be described by a power law r~ 7 / 3 in a region of size 
R s around the IMBHs. Annihilations themselves will set an upper limit to the DM density 
Pmax ~ m x/[{°~ v )Ai where t is the time elapsed since the formation of the mini-spike, and we 
denote with R c the "cut" radius where p(R c ) = /Vax- With these definitions, the intrinsic 
luminosity factor in Eq. El reads 

-^sp = PlOO (^s)^ 3 ^ (4) 

where -R SiPC and -R c ,mpc denote respectively R s in parsecs and R c in units of 10~ 3 pc, piqn£r) 
is the density in units of lOOGeV cm -3 . Typical values of L sp lie in the range 0.1 - 10 [16]. 
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The shape and normalization of the vg spectrum per annihilation N Ue (E), depends on the 
nature of the DM particle. Here we work out the prospects for detection in two different 
annihilation scenarios, one corresponding to annihilation to bb pair s, representative of a 
large class of supersymmetric models (see e.g. the discussion in |24|), and one where DM 
annihilates directly to neutrino pairs, representative of models without helicity suppression, 
such as the aforementioned in theories with Universal extra-dimensions jj- 111 the case 
of direct annihilation, we assume that all flavors are produced with equal branching ratios, 



The neutrino 



N Ue : N v : N Vt = 1 : 1 : 1, as is appropriate for Kaluza-Klein DM 
spectrum per annihilation is in this case mono-energetic N Ue (E) = 25(E — m x )/m x . The 
case of annihilation to bb is more complicated, since in this case neutrinos originate from 
the decay of hadrons produced in quark fragmentation, as well as from semi-leptonic decay 
of b quarks (see e.g. the discussion in ref. |25|). A typical fragmentation and decay chain 
is b — > 7r ± + 7r =1= — > fj, + u^, /i — > e + + v e . Here we fold pion spectra obtained from 
fragmentation functions j^J, with the spectra of neutrinos from pion decay (see e.g. j^). 
For semi-leptonic decay, sub-dominant at low values of E/m x , but inducing a hardening of 
the spectrum at high energies, we follow the treatment of Ref. 28|. 

We focus here on the flux at the ground, since neutrino telescopes are most sensitive 
to muons produced by charged- current interactions of with nuclei around the detector. 
However, we cannot use directly Eq. El because flavor ratios are modified by neutrino oscil- 
lations in vacuum. For the energies (E > lGeV) and distances (D^ pc > 1) considered here, 
neutrino oscillations will be completely averaged out 0, Q The flux at the detector can 
in this case be expressed as 

<M*o= E p("t^^X( E ) ( 5 ) 

£=e,/j,,T 

where P{yi — > v^) is the probabilit y fo r a to oscillate into v^, which depends on the 
neutrino oscillation matrix elements j^ll 

P{v^ v») = Y.\ U lW U U ( 6 ) 

j 

The analysis of solar and atmospheric neutrino oscillation, as well as a number of reactor 



and accelerator experiments, severely constrain the oscillation matrix elements 31] , and we 
adopt here P(ug — > u^) = (0.2,0.4,0.4) j^. The rate of muons induced by a neutrino flux 
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FIG. 2: Number of IMBHs producing a rate R, or larger, in Antares, IceCube and Km3. Solid 
lines corresponds to the average over all realizations, while dashed line denote the 1-a scatter from 
one simulation to another. For reference we show the rate induced by atmospheric neutrinos in 
Antares (vertical dotted) and Km3 (vertical dashed) in a 1° cone around a source located at the 
Galactic center. 

(E) is given by 

r m x rv» 

R = V 4> (5) dE u dyA(E„)P^E V} y)^{E v ) (7) 
Je^ 1 Jo 

where y v = 1 - E^ T / E v and E^ r is the muon threshold energy of the experiment (we adopt 
here _E* hr = 100 GeV). A(E lt ) is the effective area of the detector. In our case we have used 
an energy-dependent A(E^) for Antares [3^, while for IceCube we assumed the indicative 
value A(E ft ) = 1km 2 . Pp(E v ,y) is the probability that a neutrino of energy E v interacts 
with a nucleon producing a muon of energy E^ = (1 — y)E v above the detector threshold 
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energy, and can be estimated as 

P,(E U , y) = N A R(E„ E^) a(E u , y) (8) 

where R(E^ Ef T ) is the muon range, i.e. the distance traveled by muons before their 
energy drops below -E* hr , Na = 6.022 x lO 23 ^" 1 is Avogadro's number and a(E u ,y) = 
d<Jcc(E u ,y)/dy is the differential cross section for neutrino-nucleon charged-current scat- 
tering (for further details see e.g. Ref. ||| and references therein). In our treatment of 
high-energy neutrino-nucleon interactions we follow Ref. j^jj, but we make use of an up- 
dated set of parton distribution functions [35]. 

The results are shown in Fig. |21 where we plot the number of black holes producing a rate 
of events R, or larger, in ANTARES and IceCube, assuming m x = 1 TeV, and av = 10 26 cm 3 
s -1 . For comparison, we show the rate of atmospheric neutrino events in a search cone of size 
1° around the source, for Antares and for a kilometer-scale telescope. For the atmospheric 
background in Antares we have adopted the value derived in a dedicated study , relative 
to a point source coincident with the Galactic center, while for kilometer scale telescopes, 
we have used the so-called Bartol flux [3y|, and a source at the same declination. 

The choice of a heavy DM particle in Fig. El was motivated by the fact that the sensitivity 
of neutrino telescopes increases with energy, due in particular to the increase of the neutrino- 
nucleon cross-section, while the atmospheric background rapidly decreases with energy, a 
circumstance that leads to a better signal to noise ratio for high-energy fluxes. Unlike other 
indirect searches, annihilation fluxes from mini-spikes around IMBHs are not very sensitive 
to the particle physics parameters of the DM candidate. In fact, the luminosity factor L sp 
implicitly depends on the DM mass and cross-section, partially compensating the (cxt>)/m 2 
factor in Eq. El so that the neutrino flux is actually proportional to (m x )^ 9 ^ 7 (av ) 2 / 7 , as 
for gamma-rays [IfiJ. An additional reason to focus on a heavy DM candidate is that the 
prospects for detection with gamma-ray telescopes are more promising for DM particles 
lighter than 300 GeV. 

As one can see from Fig. IceCube is expected to detect a larger number of sources with 
a high event rate, with respect to Antares, and several of those are expected to exceed the 
atmospheric background. However, the sources are more concentrated toward the Galactic 
center, i.e. a region that cannot be seen by IceCube. The Antares location is more favorable, 
although the event rate is inevitably smaller due to reduced size of the experiment. A kilo- 
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meter sized detector in the Mediterranean |37|, combining the large, km 2 , area of IceCube, 
with the visibility function of Antares, appears ideally suited for this type of searches. We 
show in Fig.|21what can be achieved with a detector with A{E P ) = 1km 2 at the same latitude 
as Antares. Note also that the specific example in Fig. |2]is consistent with null searches with 
the Amanda-II telescope although the collected data can already be used to constrain 
the parameter space of the proposed scenario. 

The results for direct annihilation to V(pi are less promising, unless the mass of the DM 
particle is just above E^. This is due the fact that in the bb channel, fragmentation and 
decay of secondary particles produce a large number of neutrinos, with energy much smaller 
than m x . The direct annihilation thus dominates only when most of these neutrinos are 
produced at energies below -E* hr . For the parameters in Fig. El the vipi flux is smaller by a 
factor ~ 3.5 for Antares, and ~ 5 for IceCube and Km3. Note that although neutrinos from 
the Galactic center may in principle represent an interesting alternative, the prospects for de- 
tectog g_ays from. source are more Werestog both for superset* Q 
and Kaluza-Klein DM |25j,|4Qj. 

We recall that we derived our results in the framework of a specific IMBHs scenario, 
characterized by a heavy mass scale ~ 1O 5 M . Alternatively, we could have chosen to work 
in a different scenario, where IMBHs form in overdense regions at high redshift, z ~ 20, 
as remnants of Population III stars, and have a characteristic mass-scale of a few 10 2 

n n 

41]. However, in this scenario the formation of IMBHs is not adiabatic jlfil ]. and mini 



42, 



spikes can only form by the subsequent growth of IMBHs by accretion of baryons 
1-431 ] . a circumstance that makes the mini-spikes more shallow, and the whole scenario more 
problematic. 

We have shown that upcoming neutrino telescopes may detect the neutrino flux from DM 
annihilations around IMBHs. The prospects for detection appear promising, especially for a 
kilometer-scale neutrino telescope in the Mediterranean. Both Antares and IceCube should 
be able to detect several, up to many, IMBHs. 

Unlike searches of high energy neutrinos from the Sun or the Earth, which in most cases 
are only sensitive to scattering cross-section of DM particles off-nuclei, the proposed search 
is only sensitive to the annihilation cross section, a circumstance that makes the proposed 
search complementary to the existing ones. 

We stress the importance of a dedicated analysis that keeps into account the peculiar ex- 
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perimental details of the various neutrino telescopes, and a careful study of the detectability 
of sources in different regions of the DM parameter space. 
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